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ABSTRACT 

We study cosmic-ray acceleration in a supernova remnant (SNR) and the 
escape from it. We model nonthermal particle and photon spectra for the hidden 
SNR in the open cluster Westerlund 2, and the old-age mixed- morphology SNR 
W 28. We assume that the SNR shock propagates in a low-density cavity, which is 
created and heated through the activities of the progenitor stars and/or previous 
supernova explosions. We indicate that the diffusion coefficient for cosmic-rays 
around the SNRs is less than ~ 1% of that away from them. We compare our 
predictions with the gamma-ray spectra of molecular clouds illuminated by the 
cosmic-rays (Fermi and H.E.S.S.). We found that the spectral indices of the 
particles are ~ 2.3. This may be because the particles were accelerated at the 
end of the Sedov phase, and because energy dependent escape and propagation 
of particles did not much affect the spectrum. 

Subject headings: radiation mechanisms: non-thermal — ISM: individual (W 28) 
- cosmic rays — supernova remnants — open clusters and associations: indi- 
vidual (Westerlund 2) 



1. Introduction 



Supernova remnants (SNRs) are canonically considered the main sources of cosmic-rays 
in the Galaxy. The detection of non-thermal X-ray emissio n from SNRs clearly indicates that 



electrons are actually accelerated around the SNR shocks (IKoyama et al.lll995l ). and the ob- 
servations can constrain the electron spectra. On the other hand, observational confirmation 
of accelerated protons is not as easy as that of electrons. One way to study the acceleration 



decay of neutral pions (e.g. 


STaito & Takahara 


1994: 


Drurv et al. 


1994 




Sturner et al. 


1997; 


Gaisser. Protheroe & Stanev 


|l998; 


Baring et al 


. 19991: 


Berezhko & Volk 


2000). In particu- 



lar, molecular clouds are efficient targets of cosmic-ray protons because of their high density. 
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Thus, clouds ill uminated by the protons accelerated in a nearby SNR could be bright gamma- 



ray sources (e.g. 


Aharonian & Atovan 


1996 


Fatuzzo & Melia 


2005; 


Fatuzzo, Adams & Melia 


2006; 


Gabici. Aharonian & Casanova 


2009 


)• 



Theoretical studies have suggested that old SNRs could be appropriate objects to in- 
vestigate gamma-ray emission through pp-interactions, because the radiation from the ac- 
celerated elec trons (primary electrons) disappears as the SNR evolves, owing to their short 
cooling time (jYamazaki et al.ll2006l ; iFang fc Zhangl 120081 ) . In other words, we could ignore 
the gamma-rays from primary electrons via inverse- Compton (IC) scattering of ambient 
soft photon fields and/or non-thermal bremsstrahlung from the interaction of electrons with 
dense ambient matter. 

In this letter, we consider the evolution of an SNR surrounded by molecular clouds. We 
calculate the spectrum of cosmic-rays accelerated in the SNR and the photon spectrum of 
the molecular clouds illuminated by the cosmic-rays. We assume that a supernova explodes 
in a low-density cavity, because the progenitor star expels ambient dense gas via strong 
UV-radiation and stellar winds (jChevalierl Il999i ) . The cavity may also have been created 
through previous supernova explosions. What differentiates this study is that we consider 
whether high-energy cosmic-rays illuminating molecular clouds were accelerated even after 
the SNR became old or they were accelerated only when the SNR was young. We also discuss 
the influence of particle diffusion on the cosmic-ray spectrum. We construct specific models 
for the open cluster Westerlund 2 and the SNR W 28, and compere the results with latest 
observations. 



Westerlund 2 is one of the yo ung open clusters from which TeV gamma-ray emission 
has bee n detected with H.E. S.S. ( lAharonian et al.l 120071 ). It is surrounded by molecular 



clouds (IFurukawa et a 



2009) . The gam ma-ray emission is extended (~ 0.2°) and covers 



the molecular clouds (IFukui et al.ll2009l ). Noticeable objects such as pulsars that can be 



the so urce of the gamma-ray emission have not been observed in this region. iFujita et al. 



(120091 ) proposed that the gamma-ray emission comes from an old SNR, although there is no 
clear signature o f SNRs in the c luster. W 28 is a mixed-morphology SNR interacting with a 
molecular cloud (IWoottenlll98ll ). It is a n old SNR and TeV gam ma-rays have been detected 
from molecular clouds around the SNR (lAharonian et al.ll2008l ). 



2. Diffusion of Cosmic-Ray Particles 



As will be shown in the next section, the proton spectrum around Westerlund 2 and 
W 28 can be fitted with a power-law with an index s ~ 2.3 (the definition is shown in 
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equation [4]). Moreover, the shock waves seem to have finished accelerating particles for 
those objects, while the surrounding regions are bright in the gamma-ray band. Here, we 
briefly discuss what they mean before we explain assumptions in our specific models for 
Westerlund 2 and W 28. 



2.1. Energy Spectrum of Cosmic-Rays 

If the duration of particle acceleration is shorter than that of the di ffusion, and the parti 



cle sou rce is spatially localized well, we can use the analytical solution in lAtoyan. Aharonian. fc Volk 



( 119951 ). This corresponds to the situation where particles are accelerated mainly when the 
SNR is young and compact, and the molecular cloud illuminated by accelerated cosmic-rays 
is distant from the SNR shock. If the shape of the source spectrum is a power-law with 
an index a or Q oc E~ a 8(r)8(t), the energy spectrum at the position of the cloud (r) is 
represented by 

E~ a ( r 2 \ 

/(^)«7^ ex P -iH , (i) 

' diff V ' diflf / 

if radiative cooling during the diffusion can be ignored. The diffusion lengt h is r^s = 



2^J D{E)t ) where D{E) is the diffusion coefficient. Following iGabici et al.l f 20091 ) . we assume 
that 



where B is the magnetic field. At a point distant from an SNR, we expect that x ~ !• Thus, 
for a given magnetic field, the energy spectrum is represented by / oc E~ a ~ 75 if r < r^e 
at the position of the molecular cloud. This means that even if particles are accelerated 
efficiently (a = 2), the energy spectrum must be soft (/ oc E~ 2 ' 75 ). In other words, if the 
index of the spectrum is observed to be s < 2.75 at a molecular cloud, it is likely that the 
particles are accelerated near the cloud after the SNR becomes large. 

For Westerlund 2 and W 28, since the spectral indices are s ~ 2.3 for the high-energy 
protons (> TeV) that are illuminating the molecular clouds around these objects, we expect 
that the cosmic-rays were accelerated near the molecular clouds even after the SNRs became 
old and large. This may be in contrast with the assumption often adopted in theoretical 
studies. We assume that the SNR shock had traveled in a low-density cavity. During the 
propagation in the cavity, the shock wave is in the adiabatic Sedov phase, because the low- 
density prevents the radiative cooling. Thus, even if particles can be accelerated only during 
the Sedov phase, they are being accelerated until the shock reaches and collides with the 
surrounding high-density region, which is an effective target of the cosmic-rays. The particles 
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illuminate the high-density region with the energy spectrum at the acceleration site or the 
shock. Thus, the spectral indices of s < 2.75 are possible. 



2.2. Diffusion Coefficient 



For Westerlund 2 and W 28, efficient particle acceleration around shock waves seems to 
have finished. For the former, there is no signature of a shock. Thus, the shock may have 
already collided with the surrounding dense region and dissipated. For the latter, the shock 
is trave ling in a relatively high d ensity region (~ 10-100 cm -3 ) and seems to be in a radiative 
phase (jllho fc Borkowski 2002 ). Thus, cosmic-rays that are illuminating mol ecular clouds 



might be accelerated in the past. For x — 1 an d B = 20 fiG in equation iGabici et al. 
(120091 ) estimated the diffusion time of cosmic-rays in a typical molecular cloud and showed 
that it is only ~ 100 yr for protons with an energy of ~ 10 TeV. However, the probability 
that the particle acceleration has finished within ~ 100 yr must be very small. Therefore, 
in the next section we simply assume that x ^ 1 around the SNRs and later discuss it 
quan titatively . The small x may be because of the generation of plasma waves by cosmic- 
rays flWentzej[l974f ). 



Note that if particles are accelerated with a spectral index of a = 2, and if they diffuse 
in the interstellar gas and illuminate molecular clouds before the index increases to s — 2.75, 
the observed spectral index could be s ~ 2.3 as the ones observed in Westerlund 2 and W 28. 
However, the value of x and/or the distance between the SNR and the molecular clouds must 
be fin e-tuned in this case. In particular, for W 28, the molecular clouds are adjacent to the 
SNR (lAharonian et al.l 120081 ) and there seems to be no room to adjust the distance. Thus, 
we do not consider this possibility. 



3. Models and Results 



3.1. Westerlund 2 



We assume that the distance to Westerlund 2 is d = 5.4 kpc (IFurukawa et al.ll2009l ). We 
also assume that a supernova exploded at the center of a cavity, which is filled with warm 
gas, and is surrounded by molecular clouds (Fig. dh). After the explosion, a shock expands in 
the warm gas and then hits the surrounding molecular gas ( ~ 20 pc from t he ca vity center). 
The size of the cavity is based on the arc structure found by iFukui et al. J2009h . 



When the shock expands in the cavity, the shock velocity v s is written as a function of 
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the SNR age (t) as 



v s (t) 



Vi 

Vi{t/t X ) 



(t < h) 

" 3 / 5 (t!<t< t 2 ) 



(3) 



wher e Vj = Vj^lO 9 cm s 1 is the initial velocity of the ejecta ( Sturner et al.lll997uYamazaki et al 



20061 : iFang fc Zhand 120081 ) . At t < t x = 2.1 x 10 2 (E 51 /n 



1V3,,- 5 / 3 



J i,9 



yr, the SNR is in the 



free expansion phase, at t > t 2 = 4.0 x lO 4 ^ 



4/17^-9/17 
51 



n. 



yr, it is in the radiative phase, and 



at ti < t < t 2 , it is in the Sedov phase. Here n = 100 n 2 cm -3 and n = /i n CSLV cm -3 , where 
n cav is the hydrogen density in the cavity and /i — 1.4 is the mean atomic weight of the gas 
assuming one helium atom for every 10 hydrogen atoms. The initial energy of the ejecta is 
E SN = .E^IO 51 erg. The radius of the shock (r s ) can be obtained by integrating equation (j3J). 



The energy spectrum of the accelerated protons is given by 



N P (E) oc E~ s exp(-E/E, 



max,p J 



(4) 



The index is given by s = (r cr , m + 2)/(r com — 1), where r com is the compression ratio of 
the shock (IBlandford fc Eichlerlll987l ). which is given by the Rankine-Hugoniot relation for 
ti < t < t 2 - We do not consider non-linear effects, and thus r com is equal to or smaller than 
4. We assume that the maximum energy is determined by the age of the SNR: 



max,p 



1.6 x 10 z h 



2 u-i 



10 s cm s 



-1 



10 fiG J \ 10 5 yr 



t 



TeV 



(5) 



where h (~ 1) is the factor determined by the shock angle and the gyro-factor, and is 
the downstream m agnetic field, which is g iven by B d = r com B cav , where B cav is the magnetic 
field in the cavity (jYamazaki et al.ll2006l ). The minimum energy of the protons is given by 
the rest-mass energy. 



In our fiducial model, we assume that n c 



4 cm 3 which is the lower end of the 



typical value in such regions (~ 5-25//i cm -3 ; IChevalierl Il999f ) . There is no information 
about the temperature and magnetic field in the cavity before the explosion. Thus, we 
assume that T cav = 8 x 10 5 K and S cav = 20 jiG so that the photon spectrum is consistent 
with observations; T cav and B cav determine s and E may[tP , respectivel y. It is to be noted that 
the current temperature of the warm gas is 1.4 x 10 6 -4.9 x 10 7 K f lFuiita et all 120091 ). For 
Ui 9 = E 5 i = 1, the shock radius is r s = 20 pc at t = 1.6 x 10 4 yr. We define this time as £ co u 
and the shock collides with the molecular clouds around this time. Since t 2 = 1.6 x 10 4 yr, 
the SNR is at the end of the Sedov phase at t — t co ii. Assuming that particle acceleration 
stops at t ~ t 2 , the protons accelerated at t < t co ii (~ t 2 ) illuminate the molecular clouds. 
For ti < t < tcoii, the compression ratio r com and the proton cut-off energy £/ max ,p decrease 
from 4 to 3.6 and from 151 to 47 TeV, respectively. Therefore, the protons that illuminate 



the molecular clouds have a spectrum represented by equation with 2.0 < s < 2.3 and 



47 < £ max , p < 151 TeV. 



The mass of the molecular clouds around Westerlund 2 is ~ 2 x 1O 5 M ( IFurukawa et al. 



20091 ). although the clouds are very irregular. Referring to the observations, we assume that 
the proton density of the clouds is n c = 1.0 x 10 3 cm -3 . In the clouds, we expect that the 
shock follows the "snowplow" evolution (t > £ C oii)- We assume that the shock in the clouds 
can locally be expressed as a shell centered on r = 0. In this case, from the momentum 
conservation, the radius of the shell, r S h(t), has a relation of 



4tt 



An 



[n c (r sh (t) 3 - r s (t co n) 3 ) + n cav r s (t coll ) 3 ] r sh (t) = — n 



7 T S^ColO^S^COll) 



(6) 



with r. 



sb 



takes t stop ~ 



r s at t — t co u, which can be solved numerically for t > t co \\. After the collision, it 
1 .8 x 10 4 yr until the sh ell velocity decreases to the internal velocity of the clouds 



20 km s 1 ; IFurukawa et al.ll2009l ). Since the shock has not been observed apparently, its 



velocity may have decreased to that level. 

Some of the accelerated protons plunge into the molecular clouds and create elec- 
trons, positrons, and gamma-ray photons through pp-interactions there. The electrons and 
positrons (both are called secondary electrons) also emit photons through synchrotron, IC, 
and bremsstrahlung radiation. Figure [2] shows the broad-band spectrum of the molecular 
clouds at t ~ t co ii- We assume that the spectrum of the cosmic-ray protons is represented by 
equation PJ) with s = 2.3 and E maKtP = 47 TeV. The total energy of the protons illuminating 
the molecular clouds is 4.5 x 10 49 erg. Note that the total number of the target protons and 
the total energy of the illuminating protons are degenerated for the photon spectrum. The 
magnetic field i n the cloud is Br = 6 fJ,G. The spectrum is calculated based on the radi- 
ation models in iFang &: Zhand (120081 1 . For the pp-interactions, we used the code provided 



by Karlsson Sz Kamael ( 2008 ) . In Figure El we included the contribution of background 



cosmic-rays (equation [10] in iGabici et al.l 120091 ). but it is negligibly small. We assumed 
that the injection of secondary electrons is balanced with the cooling. The model results 
are compared with observations. The radio flux is considered as an upper limit, because it 
could only partially includ e the non-thermal radio flux produced by the energetic electrons 
(IWhiteoak fc Uchidalll997l ). For the two-band Fermi observations, we assume that the spec- 
tral index is between —2.5 and —2. The data were particularly important to constrain the 
value of s. 



The diffusion time of cosmic-rays in a cloud is written as 
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where L c is the size of a molecular cloud. From the condition of t<ns ^ t siop , the reduction 
factor must be \ < 0.01. If we ignore the diffusion of the protons in the cloud, the emission 
originated from the proto ns will not change for a long time because of the long cooling time 



of the protons (> 10 5 yr; iGabici et al.l 120091 ) . 



In Figure [21 the emission from primary electrons is also presented. Their energy spec- 
trum is given by 

N e (E) oc E~ s exp(-£/£ max , e ) . (8) 
The maximum energy is determined by synchrotron cooling: 

£max,e = 14 h^v^B^ '10 fiG)^ 2 TeV , (9) 



(jYamazaki et al.l 120061 ) . The minimum energy of the electrons is given by the rest-mass 
energy. We assume that the electron-proton ratio is K ep = 0.01. If the bremsstrahlung 
radiation from primary electrons were more than ten times brighter than that in Figure [21 
it might overwhelm the gamma-ray emission from the pion decay. However, the required 
electron-proton ratio (K ep > 0.1) would be unrealistically large. Moreover, because of the 
short cooling time, the radiation of the primary electrons should decrease after particle 
acceleration stops. Thus, for t > t co \\ (~ t 2 ), observed flux from the primary electrons should 
be lower than that in Figure [21 



3.2. W 28 



We apply the same model to W 28. A supernova explodes in a low- density cavity 
(Fig. [Tb). The distance to W 28 has been estimated to be d ~ 2-3 kpc (IGoudisI 1 19761 ; 



= 2.6 kpc in this letter. We ass ume that n nav = 8 cm 3 
K. For v ifi = 1 and E 51 = 0.4 faho fc Borkowskil 12002^ 



Lozinskayal 1 19811 ) ; we assume d - 
£ cav = 2 //G, and T cav = 1 x 10 6 
the shock radius is r s = 11 pc at t = t% = 9.0 x 10 3 yr. During the Sedov phase, the 
index s increases from 2.0 to 2.3, and the maximum proton energy -E maXjP decreases from 
8.8 to 2.7 TeV. At t ~ t 2 , we assume that the shock collides with a surrounding relatively 
dense region of rih = 70 cm " 3 (Fig. [3b) and we define this time a s t = t co a. The density 
rih is based on observations ([Long et al.lll99ll ; iReach fc Rhdl2000l ). The detection of 1720 
MH z maser emission from W 28 actually indicates the presence of the shock in the dense 
gas ([Frail et al.lll994l ; IClaussen et al.lll997l ). After the collision, the evolution of the shock 
follows equation in which n c is replaced by n^. The radius of the shell increases from 
r s h = 11 to 13 pc, and the velocity decreases from v s = 500 to 80kms _1 in t stop ~ 1.4 x 10 4 yr. 
The radius of r, s = 13 pc and the velocity of 80 km s" 1 are comparable to the observed ones 



( IBohigas et al.lll983l ; iRho fc Borkowskil 120021 ) . Note that even if particles can be accelerated 
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at radiative shocks, it cannot be applied to W 28 at present, because the small velocity 
kms -1 ) prevents ionization of the preshock gas (jShull fc McKedll979l ). 



Since the gamma-ray emission comes from clumps around the SNR (jAharonian et al. 



20081 ). we expect that there are dense clumps, which are the targets of cosmic-ray protons 



(Fig. [lb). Referring to the observations, we assume tha t the mass and density of a clump 
is 4 x 10 4 Mq and n c = 5.4 x 10 2 cm -3 , respectively (jAharonian et al.l 120081 ). Figure [3] 
shows the broad-band spectrum of the molecular cloud at t ~ t coll . The spectrum of the 
cosmic-ray protons is represented by equation (j4j) with s = 2.3 and £ maj[]P = 2.7 TeV. The 
total energy of the protons that contribute to the illumination of the molecular clouds is 
2 x 10 49 erg. The magnetic field in the cloud is B c = 60 fiG. The radio flux is considered as 
an upper limit, because it includes radio em ission outside of the region where TeV gamma- 
ray has been detected (IDubner et al.l l2000l ). Since the XMM- Newton upper-limit is for 
part of the TeV gamma-ray region, the actual upper-limit may be larger (Nakamura et 
al. 2009, in preparation). As is the case for Westerlund 2, while the photon flux from 
primary electrons decreases as the time advances, that originated from protons does not. 
The condition t<ns > ^sto P gives x ^ 0.01. Recently, Fermi fo und GeV gamma -rays in 
the molecular clouds where TeV gamma-ray has been detected (lAbdo et al.l 120091 ) , which 
indicates that the energy dependence of D(E) is not affecting the proton spectrum and is 
consistent with our assumption that the protons were accelerated near the molecular clouds. 



4. Summary and Discussion 



In this letter, we constructed a model of gamma-ray emission from molecular clouds 
surrounding an SNR. We considered old SNRs because gamma-rays mainly come through 
pp- interactions and pion decay, and those from primary electrons can be ignored. In fact, 
the synchrotron co oling; time of electrons with E ~ 1 TeV is only ~4x 10 3 yr for B ~ 60 fiG 
(IGabici et al.ll2009l ). and thus they cool quickly after their acceleration finishes. We assumed 
that the SNR shock travels in a low-density cavity, and we indicated that the diffusion 
coefficient of cosmic-ray particles around the SNR needs to be much smaller than that away 
from it (less than ~ 1%), The predicted broad-band spectra for Westerlund 2 and W 28 
are consistent with observations (Figs. [2] and [3]). The proton spectra can be fitted with a 
power law with indices ~ 2.3, which i s the value at cosm ic-ray sources expected in cosmic-ray 
propagation models for the Galaxy (jStrong et al.l 120071 ) . The indices suggest that particles 
were accelerated near the molecular clouds even after the SNR became old and large for 
these objects. The soft spectra may be because the particles were accelerated at the end 
of the Sedov phase. While the low Mach number of the shock makes the spectra soft in 
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our m odel, neutral hydrogens in the cavity may also do that (jOhira. Terasawa. fe Takahara 



20091 ). On the other hand, for another SNR IC 443, the broad-band photon spectrum 



indicates that the index is 2 (jZhang fe Fang||2008l ). For this object, protons accelerated at 



an earlier stage of the Sedov phase (when the Mach number of the shock is large) may be 
illuminating molecular clouds. Moreover, if molecular clouds are distributed with various 
sizes and distances from the SNR, the gamma-ray spectrum may be the superposition of 
the spectrum of each cloud and may reflect the energy-depe ndence of particle diffusion 



(equation [2J; see lAharonian fe Atoyanlll996l ; iTorres et al.ll2008l ). 
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Fig. 1. — Configuration of the cavity and clouds, (a) Westerlund 2 and (b) W 28. 



13 



10" 



'co 10 -10. 




Fermi 
Suzaku jc^l ■ 



10 L 
E(MeV) 



w 
w 

\\\ 

U J 



10 



8 



Fig. 2. — Compar ison of t he model results w ith radio (IWhiteoak fc Uchidalll997T). Suzaku 



fFui ita et al.ll2009l ). Fermi (lAbdo et al.ll2009l ). and H.E.S.S. observations (lAharonian et al. 
2007) for Westerlund 2. The synchrotron radiation (solid line), bremsstrahlung (dashed line), 
and IC scattering (dotted line) are of the primary electrons (thin lines) and the secondary 
electrons (thick lines). The 7r° decay gamma-rays are shown by the two-dot-dashed line. 




Fig. 3. — Comparison of the model results with radio (IDubner et al. 120001). XMM- New ton 
(Nakamura et al. 2009, in preparation), and Fermi observations (lAbdo et al.l 120091 ) for 



W 28. H.E.S.S. observations are for HESS J1801-233 jAharonian et al.lboosl ). The lines are 
the same as those in Fig. [2j 



